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We discuss the physics potential of intermediate L ~ 20 30 km baseline experiments at reactor 
facilities, assuming that the solar neutrino oscillation parameters ArriQ and 6q lie in the high-LMA 
solution region. We show that such an intermediate baseline reactor experiment can determine 
both AmQ and 9q with a remarkably high precision. We perform also a detailed study of the 
sensitivity of the indicated experiment to Arriatm, which drives the dominant atmospheric {p^) 
oscillations, and to 6 - the neutrino mixing angle limited by the data from the CHOOZ and Palo 
Verde experiments. We find that this experiment can improve the bounds on sin^ 9. If the value 
of sin^ 9 is large enough, sin^ 9 ^ 0.02, the energy resolution of the detector is sufficiently good 
and if the statistics is relatively high, it can determine with extremely high precision the value of 
Am^tnj. We also explore the potential of the intermediate baseline reactor neutrino experiment 
for determining the type of the neutrino mass spectrum, which can be with normal or inverted 
hierarchy. We show that the conditions under which the type of neutrino mass hierarchy can be 
determined are quite challenging, but are within the reach of the experiment under discussion. 

PACS numbers: 14.60.Pq 13.15.+g 



I. INTRODUCTION 

The experiments with solar, atmospheric and reactor neutrinos Q, |M 1^ Q IM IS have provided in recent years 
remarkable evidences for the existence of neutrino oscillations driven by nonzero neutrino masses and neutrino mixing. 
The hypothesis of solar neutrino oscillations, which in one variety or another were considered as the most natural 
explanation of the solar neutrino deficit 0, Q since the late 60'ies (see, e.g., 0, H 0), has received a convincing 
confirmation from the measurement of the solar neutrino flux through the neutral current reaction on deuterium by 
the SNO experiment 0- The analysis of the solar neutrino data obtained by Homestake, SAGE, GALLEX/GNO, 
SK and SNO experiments showed that the data favor the Large Mixing Angle (LMA) MSW solution with the two- 
neutrino oscillationparameters - the solar neutrino mixing angle and the mass squared difference, lying at 99.73% 
C.L. in the region |alTol|: 

3 X lO^^eV^ < Am^j < 3.5 x lO^^eV^ (1) 
0.21 < nin^OQ < 0.47. (2) 

The first results of the KamLAND reactor experiment has confirmed, under the plausible assumption of CPT- 
invariance which we will suppose to hold throughout this study, the LMA MSW solution, establishing it essentially 
as a unique solution of the solar neutrino problem. The combined fits to the available solar neutrino and KamLAND 
data, performed by several collaborations within the two-neutrino mixing hypothesis, identify two distinct solution 
sub- regions within the LMA solution region 0, IT^ [Ts j . Adding the KamLAND data did not lead to a considerable 
reduction of the interval of allowed values of sin^ 9q with respect to the one quoted in eq. 0J , while the best fit values 
of Atoq in the two sub-regions (labeled from now on low-LMA and high-LMA) are given by '■ 

low-LMA: Am| ~ 7.2xlO-^eV^ (3) 
high-LMA: Am| ~ l.SxlO^^eV^ (4) 

The observed Zenith angle dependence of the multi-GeV /i— like events in the Super-Kamiokande experiment un- 
ambiguously demonstrated the disappearance of the atmospheric (D^) on distances L ^ 1000 km. The Super- 
Kamiokande (SK) atmospheric neutrino data is best described, as is well-known, in terms of dominant j/^ (y^^ — > 



Vr) oscillations withf almost) maximal mixing and neutrino mass squared difference of |Am^j,j^| = (1.4 — 5.0) x 10 
eV^ (99.73% C.L.) According to the more recent combined analysis of the data from the SK and K2K experiments 
fil one has at 99.73% C.L.: 



1.5 X 10"^ eV^ < IAt 



'atml 



< 3.7 X 10"^ eV^ 



(5) 



The neutrino oscillation description of the combined solar and atmospheric neutrino data requires the existence of 
three-neutrino mixing in the weak charged lepton current (see, e.g., {15]): 
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Here i/j^, I = e,iJ,,T, are the three left-handed flavor neutrino fields, z/^l is the left-handed field of the neutrino Vj 
having a mass rrij and U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mixing matrix [ll[l3. The 
PMNS mixing matrix U can be parameterized by three angles, ^atm, ^o, and 9, and, de pending on whether the 
massive neutrinos Vj are Dirac or Majorana particles - by one or three CP-violating phases [iMIl^- In the standard 
parameterization of U (see, e.g., 15j) the three mixing angles are denoted as O12, O13 and ^23: 
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where we have used the usual notations, Sij = sin0. 



, = cos 6ij , and 5 is the Dirac CP- violation phase ^ . If we 
identify the two independent neutrino mass squared differences in this case, Am|]^ and Am§]^, with the neutrino mass 
squared differences which drive the solar and atmospheric neutrino oscillations, Amg — Am^i > 0, Am^j,-,-^ — Arri^i, 
one has: 6*12 = 9q, 623 = 6'atm, and 6*13 = 9. 

The angle 9 is limited by the data from the CHOOZ and Palo Verde experiments [2lll23| which searched for evidences 
for oscillations of reactor at ~ 1 km from the source. No disappearance of Pe was observed. In a two-neutrino 
oscillation analysis performed in |21j a stringent upper bound on the value of 9 in the region of |Ato'^| > 1.5 x lO'^^eV^ 
was obtained. A S-iy oscillation analysis of the CHOOZ data ^ 23] led to the conclusion that for AttIq ^ 10^"* eV^, 
the limits on sin^ 9 practically coincide with those derived in the 2-v oscillation analysis in 21]. A combined 3-z/ 
oscillation analysis of the solar neutrino, CHOOZ and the KamLAND data, performed under the assumption of 
Am% < JAm^^^J (see, e.g., [Hi^), showed that [ij 



sin^ 9 < 0.05, 99.73% C.L. 



(8) 



The precise upper limit in Eq. (jHl is Am^^jj^-dependent. The authors of Tl'] found the best-fit value of sin^ 9 to lie in 
the interval ^ sin^ 9 ^ (0.00 - 0.01). 

Somewhat better limits on sin^ 6* than the existing one can be obtained in the MINOS experiment |2^. Various 
options are being currently discussed (experiments with off-axis neutrino beams, more precise reactor antineutrino 
and long baseline experiments, etc., see, e.g., ^27] ) of how to improve by at least an order of magnitude, i.e., to values 
of ~ 0.005 or smaller, the sensitivity to sin^ 9. 

Let us note that the atmospheric neutrino and K2K data do not allow one to determine the sign of Am^tm- This 
implies that if we identify Am'^^^ with Arri^i in the case of 3-neutrino mixing, one can have Am^i > or Am^i < 0. 
The two possibilities correspond to two different types of neutrino mass spectrum: with normal hierarchy (NH), 
nil < 1^2 < and with inverted hierarchy (IH), 1113 < mi < TO2. 

After the spectacular experimental progress made in the last two years or so in the studies of neutrino oscillations, 
further understanding, in particular, of the structure of the neutrino masses and mixing, of their origins and of 
the status of the CP-symmetry in the lepton sector requires a large and challenging program of high precision 
measurements to be pursued in neutrino physics. One of the first goals of this program is to improve the precision 



^ We have not written explicitly the two po ssible Majorana CP- violation phases llSL llfll which do not enter into the expressions for the 

oscillation probabilities of interest ll8ll2Cl . We assume throughout this study < 612, 023, S13 < 7r/2. 
^ In this case Am^ = Am^^^^^. 

^ The possibility of large sin'^ 9 > 0.97, which is a dmi tted by the CHOOZ data alone, is incompatible with the neutrino oscillation 
interpretation of the solar neutrino data (see, e.g., |25|). 



in the measurement of the mass squared differences and mixing angles which control the solar and the dominant 
atmospheric neutrino oscillations: Atoq, 0q, and Arn^^jj^, Oatm- A step of fundamental importance would be the 
detection and the studies of sub-leading neutrino oscillation effects, if the latter are observable. This includes the 
measurement of, or getting more stringent upper limit on, the value of the third and the only small mixing angle 
in the PMNS matrix U, 9 (= 6*13), the exploration of the possible CP-violating effects and the determination of the 
type of the neutrino mass spectrum which can be with normal hierarchy (NH) or inverted hierarchy (IH). Among the 
further fundamental open questions, which cannot be answered by studying neutrino oscillations, but the progress in 
the studies of which requires a precise knowledge of the neutrino oscillation parameters, are i) the nature — Dirac or 
Majorana — of the massive neutrinos, ii) the absolute scale of neutrino masses, iii) the mechanism giving rise to the 
neutrino masses and mixing, iv) the possible relation between CP-violation in the lepton sector at low energies and 
the generation of the baryon asymmetry of the Universe by the leptogenesis mechanism, v) lepton number violation 
and its possible manifestation in charged lepton flavor violating decays, to name a few. 

In the present article we explore the possibility of performing a high precision measurement of some of the 3- 
neutrino oscillation parameters in an intermediate baseline, L ~ 20 ^ 30 km, reactor neutrino experiment. We 
first address the issue of precision determination of the solar neutrino oscillation parameters, Attiq and 9q. The 
KamLAND experiment, as discussed by several authors [23, has a remarkable sensitivity to Atoq in the low-LMA 
region. However, the sensitivity of KamLAND to the value of 9q is not found to be equally good. Even under the most 
optimistic conditions KamLAND is not expected to substantially reduce the region of allowed values of 9q obtained 
from the analysis of the solar neutrino data |2fl| | . The best conditions for precise determination of 9q is in a reactor 
experimental set-up, where the baseline is tuned to a Survival Probability MINimum (SPMIN) ^. With a baseline of 
about L ^ 160 km and Am| ~ 7.2 x 10~^ eV^ or 1.5 x 10~^ eV^, KamLAND is essentially sensitive to a Survival 
Probability MAXimum (SPMAX). The sensitivity to SPMAX gives KamLAND the abihty to determine Am| with a 
high precision through the measurement of the distortion of the final state e"*"— spectrum, provided Atoq < 2 x 10^"* 
eV^. At the same time it reduces KamLAND's sensitivity to the solar neutrino mixing angle ^29]. As was shown in 
|29| |. a 70 km baseline reactor experiment could determine the solar neutrino mixing angle sin^ 9q (tan^ 9q) to within 
9.6 (14)% at the 99% C.L. in the case of the low-LMA solution of the solar neutrino problem. KamLAND essentially 
looses sensitivity to AttIq for Amg ^ 2 x 10"'' eV^. Nevertheless, the more precise measurement of the spectrum 
of the final state e+ in the KamLAND experiment is expected to unambiguously determine whether Attiq lies in the 
low-LMA or high-LMA region. 

In this paper, following a previous suggestion and analysis in [^^l (see also |23)> the context of three neutrino 
mixing, we discuss the physics potential of intermediate 1/ ~ 20 -^ 30 km baseline experiments at reactor facilities. We 
show that such an intermediate baseline reactor experiment can determine both Atoq and 9q with a remarkably high 
precision ij the solution of the solar neutrino problem is the high-LMA solution. We perform also a detailed study of 
the sensitivity of the indicated experiment to the parameters Am^tm and sin^ 9. We show that if the energy resolution 
of the detector is sufficiently good and if the statistics is relatively high, one can choose small enough energy bins so 
that even an experiment with an intermediate baseline of 20 — 30 km can be sensitive to the Am^^j^ driven oscillations. 
We find that this experiment can certainly improve the bounds on sin^ 9 and if the value of sin^ 9 is large enough, it 
can determine with extremely high precision the value of Am^^j^. 

We finally explore the potential of the intermediate baseline reactor neutrino experiment for determining the type 
of the neutrino mass spectrum which can be with normal or inverted hierarchy. The knowledge of the neutrino mass 
hierarchy is of crucial importance, in particular, for modeling of the neutrino mass matrix and understanding of the 
underlying physics of neutrino mass generation. The type of the neutrino mass hierarchy has been also shown to be 
an important "parameter" in a number of neutrino mixing and neutrino oscillation observables, such as the effective 
Majorana mass in neutrinoless double-/? decay j^^l- The neutrino mass hierarchy (or, e.g., the si gn o f Am^^^) can be 
determined in very long basehne neutrino oscillation experiments at neutrino factories (see, e.g., [23)) a-^id, e.g, u sing 
combined data from the long baseline oscillation experiments JHF-SK and NuMI with off-axis neutrino beams |34j . 
It was suggested in that the "interference" effects between the Attiq and Ato^^j^ driven oscillations can be used 
in reactor experiments to answer the question about the neutrino mass hierarchy. We show that the conditions under 
which the type of neutrino mass hierarchy can be determined are rather challenging, but are are not out of reach 
of the experiment under discussion, provided the solar neutrino oscillation parameters lie in the high-LMA solution 
region. 



* A detailed discussion of SPMAX and SPMIN as well as 60 sensitivity of the current and future experiments can be found in l2Hl and 
we do not repeat it here. 



II. THE THREE-GENERATION SURVIVAL PROBABILITY 



The expression for the i^e survival probabiUty in the case of 3 flavor neutrino mixing and neutrino mass spectrum 
with normal hierarchy (NH) is given by ^ [l^, 113, 113 : 

= 1-2 sin^ d coi? 61 I 1 - cos .^"^atm ^ 1 



IE, 

1 . , / A™?, L 



- cos^ Q sin' 29^ \ l - cos j (9) 
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where E^, is the energy. If the neutrino mass spectrum is with inverted hierarchy (IH), the Ve survival probability 
can be written in the form 12^^ 
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The survival probability does not depend either on the angle 9atm associated with the atmospheric neutrino 
oscillations, nor on the CP violating phase S in the PMNS matrix. 

In the convention (we will call A) in which the neutrino masses are not ordered in magnitude and the NH neutrino 
mass spectrum corresponds to mi < m2 < ms, while the IH spectrum is associated with the ordering < mi < 1712, 
it is natural to choose 

Am| = Am2i > , Convention A. (11) 

We can identify further Arri^^^ with Arri^i in the case of NH spectrum, 

Aml^^ = Am-li > , NH spectrum (A) , (12) 

and with Arn^^ > if the spectrum is of the IH type, 

Am^tm = ^^23 > : IH spectrum (A) . (13) 

In this convention the mixing angles in the standard parameterization of the PMNS mixing matrix U are given by: 

O12 = Oq, 023 = 6latni' O13 ^ & , Convention A . (14) 

One can also number (without loss of generality) the neutrinos with definite mass in vacuum Vj in such a way that 
their masses obey mi < 1712 < m^. In this alternative convention (we will denote as B) it is convenient to choose 

^™atm ~ ^™3i ' Convention B . (15) 



The Earth matter effects are negUgible for the values of the neutrino oscillation parameters (Arrig and Am^jj^^), Ue energies and the 
short baseline L ~ 20 30 km we are interested in. 
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FIG. 1: The two possible neutrino mass spectra, with Normal Hierarchy (NH) and Inverted Hierarchy (IH), in the convention 
in which mi < m2 < ma (convention B). 



In the case of NH neutrino mass spectrum we have 

Ato| = Am2i > , NH spectrum (B) , (16) 

and 

ei2 - 023 = 0atm> 013 = e , NH spectrum (B) , (17) 

where 9ij are the angles in the standard parameterization of U , Eq. (|7|l. 

If the neutrino mass spectrum is with IH, one has in this convention (see, e.g., |33.l36|): 

Am| = Am32 > , IH spectrum (B) , (18) 

and 



I C/e2 1 = cos Oq v/1- |C/cip, I f/c3 1 = sin v/1- iC^cip , IH spectrum (B) . (19) 

The mixing matrix element of U constrained by the CHOOZ and Palo Verde data is now |C/ciP : 

|C/eip = sin^ e , IH spectrum (B) . (20) 

We would like to emphasize that Eqs. 10 and H1U|) are valid in the two conventions for the ordering of the neutrino 
masses discussed above. 

Few comments are in order. The v^. survival probability depends only on the four continuous parameters Atoq, 
sin^^Q, Am^jjjj, sin^ 0, and on a single "discrete" parameter — the type of the neutrino mass spectrum — NH or 
IH. The terms in the second lines of Eq. ((HJ and Eq. Hl()|l are responsible for the dominant effects in the survival 
probability for the intermediate baseline experiment under consideration in the case of the high-LMA solution of 
the solar neutrino problem. The study of the dominant oscillations can thus constrain the solar neutrino oscillation 
parameters. In a complete 3-neutrino oscillation analysis, one should also consider the effect of the terms in the first 
and third lines of Eq. and Eq. (|10(l . While the terms in the first lines give the sub-dominant, faster oscillations 
driven by Am^^^^^, the terms in the last lines are responsible for the difference between the probabilities corresponding 
to the NH and IH neutrino mass spectra. Both of them are suppressed by sin^ 9. The detection of this sub-leading 
oscillatory behavior, controlled by the atmospheric neutrino mass squared difference, would be an indication of the 
non-vanishing of the mixing angle 9 in the PMNS matrix U . It will also enable the experiment under discussion to 
determine Am^^jj^ with a high precision. Finally, for sufficiently large values of sin^ 9 it should in principle be possible 
to distinguish the normal from the inverted hierarchy spectrum. Let us stress that the difference of the survival 
probabilities for the NH and the IH spectra lies in the interference term: the existing solar neutrino and KamLAND 
data indicate that 6*© is not maximal, so that sinfl© ^ cos9q. In ideal conditions in which Am^tj^ is known to a very 
high precision, sin^ 9 is not too small, sin^ 29q ^ 0.90, and the statistics and the energy resolution of the experimental 
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1.0 
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TABLE I: Number of events in the case of absence of i/^ oscillations for the different reactor experimental set-ups considered 
in this paper. 



apparatus are such to permit to resolve the Am^^j^^— driven oscillations, this could allow one to extract information 
about which of the two hierarchical neutrino mass patterns is realized in nature. 

One can rewrite the term in parenthesis in the third line of Eq. JHJ and Eq. in the form: 

^-^V-s^^ . 2 sin ^ sin f - ^) (21) 

Thus, it consists of an oscillating function, with approximately the same frequency as the Am^^.^— driven oscillations, 
modulated with period half of the usual Atoq— driven oscillations. In particular, the amplitude of this term is the 
maximal possible at values of L/Ey where the survival probability goes to its Amg— induced minima, and vanishes 
at its local maxima. 




III. THE EXPERIMENT 



Anti-neutrinos from nuclear reactor sources are detected through their inverse /3-decay reaction with protons in the 
detector. The visible energy i?vis of the emitted positron is related to the energy of the incoming anti-neutrino E^, 
and to the masses of the proton, neutron and positron as: 

E^-^^ = El, + TOe+ - (m„ - trip), 

El, - 0.8 MeV. (22) 

The i^vis spectrum of detected events in the no oscillation case, which takes into account the initial spectrum of 
anti-neutrinos emitted by the reactor and the inverse /3-decay cross section, has a bell-shaped distribution, centered 
at about E'vis ^ 2.8 MeV. The total number of expected events for the no oscillation case is given in Table We 
consider a single reactor plant with power of either 5 GW (achievable, e.g., at Heilbronn, Germany jsj) or 25 GW 
(achievable at Kashiwazaki, Japan), as the only relevant source of neutrinos, neglecting the possible contaminations 
due to other plants at larger distances. We also assume that the reactors have a 100% efficiency. It is trivial to adapt 
our results to a given reactor efficiency. We present our results as function of the product of the exposure time and 
the active mass of the detector. We will assume for this product values in the range 3 -;- 5 kTy and consider baselines 
20 30 km. The total statistics depends on C = P M T, the product of the reactor power (P), the detector mass 
(M) and time of exposure (T) , thus we express our results in units of GWkTy. Since the neutrino flux decreases as 
the inverse square of the baseline length, the shorter baselines obviously have much more statistics than the longer 
ones. We assume a liquid scintillation detector similar to ones used in the other reactor experiments like CHOOZ 
and KamLAND. We assume 8.48 x 10'^^ free protons per kton of detector mass as in KamLAND 0. We use an 
energy resolution of a{E)/E = 5%/Ve, E in MeV for "our" detector, assuming some improvement with respect to 
KamLAND, which reported a <t{E)/E ^ 7.5%/Ve 6]. 

In TableQlwe present the number of no oscillation events for two different visible energy thresholds of 1.0 MeV and 
2.6 MeV. While at higher energies the spectrum is known with relatively high accuracy, at lower energies a possible 
contamination from geophysical neutrinos and from the time variation of fuel composition is expected |37l |. The 
KamLAND experiment puts a conservative lower threshold of -Eyis '^■^ MeV to avoid the error associated with 
the geophysical neutrinos. However, since the reactor flux is the highest around 2.8 MeV, it is desirable to include 



these parts of the energy spectrum to increase the statistical power of the experiment. We note from Tabled that 
the number of observed events goes up by a factor of 1.4 when the effective threshold energy Eth is lowered from 2.6 
MeV to 1.0 MeV. Further, with a larger energy window to be used in the data analysis, and correspondingly a large 
interval in L/E^, it is possible, in principle, to reconstruct a larger number of Am^^jj^-drivcn oscillation periods, and 
even to detect both the SPMIN and SPMAX connected to the Amg-driven oscillations. The latter is crucial for the 
study of sub-leading effects performed in the second part of this paper. 

The ratio of the number of events expected from geophysics sources to the total number of antineutrinos detected 
can be estimated for a given experimental set-up. In their first published data @ , KamLAND estimated the number 
of these spurious events to be 9, while the total data set contained 32 events with energy below the energy cut-off 
{Eth = 2.6 MeV) and 54 above this threshold. For KamLAND this ratio therefore comes up to be about 30% below 
2.6 MeV. However, the number of background events due to the geophysical neutrinos is quite model dependent and 
a very good understanding and modehng (see, e.g., js^) of the energy spectrum of these events would allow one to 
use the whole data sample in the range (1.0 — 7.2) MeV, by subtracting the geophysics background from the observed 
event spectrum, with just an additional source of systematics in the analysis - the extent of the uncertainty depending 
on the accuracy of the knowledge of the geophysical neutrino background. 

The whole e"*"— energy spectrum could also be used in a very high statistics experiment with a sufficiently powerful 
reactor source. The contamination of the data sample from geophysical neutrino background at low energies is 
proportional to the active mass of the detector. Thus, decrease in the detector volume compensated with an increase 
in the reactor antineutrino flux, results in a reduction of the fraction of events due to this background. We can 
determine the flux of from a reactor in the absence of oscillations as a function of the distance traveled and of the 
power of the source: 

(23) 

For KamLAND, summing over all the reactors which contribute (see for instance for a complete listing), one gets 

^Kam _ 3 ^ Gw/ki^2 ^ ^^4) 

As a comparison, using three possible choices for the pair (L/km, P/GW) as (30, 5), (20, 5) and (20, 24), one gets for 
$ values that are bigger by a factor of 1.5, 4 and 20, respectively. Accordingly, the fraction of background events 
from the geophysical neutrinos can be reduced by this factor, as a consequence of the larger flux. In particular, for 
the last case, the geophysical neutrinos contribution would be at the percent level even for a 1 kiloton detector like 
KamLAND and could be safely accounted for. For lower fluxes and large detectors, either a very accurate subtraction 
procedure needs to be applied, or it would be necessary to keep the cutoff in the low energy part of the spectrum as 
in KamLAND. In this paper we neglect throughout the contribution coming from geophysical neutrinos - we either 
put an energy threshold of 2.6 MeV, or when we lower Eth to 1.0 MeV we assume that either they can be accounted 
for, or they can be safely neglected. 

Backgrounds may also result from radioactive impurities and from cosmic ray interactions. However, these can be 
effectively suppressed in realistic experiments and we completely neglect backgrounds in this paper. 

In presence of neutrino oscillation, the detected energy spectrum at a distance L from the reactor is obtained by 
convoluting the survival probability given by Eq. ij^ or Eq. (|10|) with the spectrum obtained with no oscillation. This 
leads to a suppression of the number of detected events, and to a distortion of the energy spectrum itself, whose 
dependence on the parameters entering the survival probabilities may be used to extract them through a fit of the 
data. 

The optimal distance in order to resolve the first minimum (SPMIN) of the survival probability, and extract with 
good accuracy a measurement of Attiq and sin^ 9q , according to Eq. Q and Eq. Hl()|l , is given by ^ : 

so that its choice depends on the actual value of AttIq and on the typical energy involved. If E^, is at the peak of the 
spectrum (3.5-^4 MeV), this choice implies a maximal depletion of the total number of events. We will consider always 
the case in which a relatively large statistics is accumulated, so that the main source of information is the differential 



This distance is also the optimal one for enhancing the effect of the interference term in the survival probability distinguishing the NH 
and IH spectrum cases, as seen from Eq. 1211 . because for this values of L/Ei, the modulation of the beating-like term is maximal. 




FIG. 2: The observed e+— energy spectrum (in 0.1 MeV bins) for the 20 km (right-panel) and 30 km (left-panel) experiment with 
reactor power 5 GW and exposure of 3 kTy (15 GWkTy). The dotted black lines in the top panels give the number of events 
in the absence of oscillations. The solid lines give the events for i?e— oscillations with ArriQ = 1.5 x 10~* eV^, siii^ 9q — 0.3, 
Am^t^ = 2.5 X 10"^ eV^. The solid cyan/ grey line corresponds to the spectrum for sin 6 = 0.0, while the solid black line 
shows the spectrum for sin^ 9 — 0.03. The bottom panels show the corresponding ratio of events in the cases of i/e— oscillations 
and of absence of oscillations for the two values of sin^ 9. The event spectra in the case of —oscillations are for neutrino mass 
spectrum with normal hierarchy. 



energy spectrum of detected events. In this case a large window of different energies can be used to determine L*. 
However, as long as the baseline of the experiment is large enough to include the first oscillation minimum (SPMIN) 
in the interval of energies used for the measurements, preferably in the part of the energy spectrum with highest 
statistics, the experiment is expected to display a good sensitivity to the solar neutrino oscillation parameters. 

We show in Fig. [3 the spectrum of observed events for two specific baselines of 20 km and 30 km for reactor power 
of 5 GW and 3 kTy exposure. We note that the 30 km experiment has much less events than the 20 km experiment. 
The solid lines in the bottom panels show that the SPMIN for the 30 km experiment comes around 2.8 MeV, where 
the flux is the highest. For the 20 km experiment the SPMIN comes at around 1.6 MeV. Correspondingly, if the 
lower energy cut-off is set to Eth = 2.6 MeV (shown by the dashed lines in Fig. 2), the 20 km experiment would miss 
the SPMIN. Therefore for the 20 km experiment, the Eth has to be lowered to below 1.6 MeV if the value of AttIq 
happens to be 1.5 x 10~^ eV^. For higher values of Attiq, the SPMIN will shift to higher energies and the threshold 
could be taken higher: for Ato| = 2.5 x lO'^ eV^, for instance, SPMIN is at E^is = 3.2 (5.3) MeV for L = 20 (30) 
km. 

Figure 121 also shows the "sub-dominant" oscillation effects dependent on Ato^^j^ and sin^ 9, imprinted on the large 
dominant oscillation wave driven by the solar neutrino oscillation parameters. The experiment could be expected to 
have sensitivity to the parameters Am^^^ and sin^ 9, driving the sub-dominant oscillations, if the energy resolution of 
the detector would be good enough to enable a large statistics experiment like the one we consider here, to break the 
total number of events into fine energy bins. In particular, to resolve the oscillatory behavior one needs bins much 
smaller than the oscillation period 

47r£'2 

AE « AE* = -— ^ , (26) 
Am^ L 

where Am^ is the relevant neutrino mass squared difference driving the oscillations. This implies that for a given 
L, higher energy resolution and smaller bins are required in order to resolve effects due to a larger mass squared 
difference. As a numerical example. 

Ami = 1.5 X 10"'* eV^ AE* = 7.1 MeV, 



Aml^ = 2.7 X 10-3 eV^ =^ Ai?*^^ = 0.4 MeV, 

for energies at the peak of the spectrum {E^ ~ 3.6 MeV) and at the distance i = 30 km. It is clear that energy 
bins with a width of 0.425 MeV (as in KamLAND) aUow a good reconstruction of the spectral distortion due to 
Attiq— driven oscillations, but would average out completely the sub-leading Ato^j^ effects. Using as a guide rule the 
assumption that to resolve an oscillation one needs Ai?^jg ^ AE* /i, the reconstruction of Amatm^driven oscillations 
requires to use energy bins having width of 0.1 MeV . The use of shorter distances can soften this problem, as 
Eq. H26() indicates. However, with an energy resolution of a{E)/E = 5%/Ve, it should be possible to have bins of 
size of 0.1 MeV, at least as long as the energy is not very large. We will show that with bins of 0.1 MeV width, and 
L = 20 km, one can improve on the CHOOZ limit for sin^ 9, and if sin^ 9 is found to be large enough to be detectable 
in this experiment, Arri^^^ can be determined with a very high precision. 

With an experimental set-up which allows a sufficiently small bin size and large enough statistics to measure the sub- 
dominant oscillations, one can hope to gain information about the neutrino mass hierarchy using the last interference 
terms in Eq. © and Eq. l(Tn|l . 

Finally, one has to take into account the presence of systematic errors both in the determination of the flux 
normalization and in the energy calibration of the detector. The KamLAND collaboration estimated a total systematic 
uncertainty of about 6.42% in their first published analysis @. The bulk of this error comes from uncertainty in the 
overall flux normalization. The use of a near detector can improve this flux normalization error to as low as 0.8% 
[3l| . The error on energy calibration could be taken as 0.5% We use a total systematic uncertainty of 2%, and 
we discuss in detail the effect of varying this uncertainty in the analysis which follows. 



IV. PRECISION NEUTRINO OSCILLATION PHYSICS : Am| AND sin^ 26*0 IN THE HIGH-LMA 

REGION 

In this Section we study quantitatively the precision measurements of the solar neutrino oscillation parameters 
Attiq and sin^ 9q , possible in the intermediate baseline reactor experiment discussed in the previous section under 
the assumption that Amg lies in the high-LMA region. We simulate the "data" at certain plausible values of Attiq 
and sin^ 9q and obtain allowed regions in the parameter space from a dedicated analysis. For the errors we assume 
a Gaussian distribution and define our as 



i,3 



where {a = data, theory) is the number of events in the z**^ bin, afj is the covariant error matrix containing the 
statistical and systematic errors and the sum is over all bins. We compare the sensitivities obtained for two possible 
baselines, 20 km and 30 km. To quantify the sensitivity we define the relative precision pa for a certain parameter a 
at a given C.L. as 

Pa = , (28) 

^max I a^nin 

where amax (omm) are the maximal (minimal) allowed value of a found at the chosen C.L. We check quantitatively 
the impact of (1) the energy threshold Eth, (2) bin size AE and (3) systematic uncertainties. We also show the impact 
of increasing the statistics by a factor of 5. 

As has been noted earlier in the context of KamLAND ||4i| , the sensitivity to sin^ 9q can be reduced considerably 
by the uncertainty in the parameter sin^ 9. We therefore analyze first the impact this uncertainty can have on the 
precision of sin^ 9q determination in the intermediate baseline reactor experiment of interest. Since the baseline is 
optimized to measure predominantly the Amg— driven oscillations, we assume that the Am^^j^j— driven oscillations 
are averaged out. This corresponds to the realistic case of using a sufficiently large e"*"— energy bin size. The expression 
for the survival probability would then reduce to 

/ Am^ L\ 

Pee W COS^ 9[l~ SIT? 29q siu^ ^® j , (29) 



where we have neglected the term ^ sin^ 9. Therefore the uncertainty in sin^ 9 essentially brings up to a ~ 10% 
uncertainty in the survival probability. Since the factor cos^ 9 can only reduce the survival probability, it does 
not affect the upper limit on the allowed range of sin^^g. However, it can have an effect on the lower limit on 



sin^ 9q reducing it further, and thus can worsen, in principle, the precision of the experiment. Using eq. H29|l we get 
approximately for the additional error on sin^ 26q due to the uncertainty in the value of sin^ 



2AP,eSin^g , . (l-sin^2go sin^ ^) A(sin^0) 
Sin' sm2 

where APee and A(sin^ 9) are the uncertainties in the determination of the survival probability and sin^ 6, respectively. 
We work in a scenario in which one would use the SPMIN in order to determine sin^ 26 q . In the SPMIN region one 
has sin^(ATOQi/4£'^) 1 and therefore 

5{siT? 29q) « 2APee sin^ 61 + 2 cos^ 26*0 A(sin2 9) . (31) 

Thus, in this SPMIN scenario, the first term gives an extra contribution of about 2APee sin^ 9 to the allowed range 
o/sin^ 26*0. Since, as we will see in the present Section, in the experiment under consideration the allowed range of 
APee ^0.1 even under the most conservative conditions, this term gives an extra contribution to the allowed range 
which is ^ 0.01. The second term is independent of the precision of a given experiment (i.e., on APee) and depends 
only on the best-fit value of cos^ 29q and on the error in sin^ 9. For the current Scr error in sin^ 9 of 0.05 and best-fit 
cos^ 29q of 0.16, this would give an increase of only 0.018. The suppression of this term is mainly due to the presence 
of the factor cos^ 26'© which is a relatively small number for the current best-fit solution. Thus, even though the 
uncertainty in sin^ 9 brings a 10% uncertainty in the value of Pee, it increases the allowed range of sin^ 29q only by 
^ (2 — 3)%, if one uses the SPMIN region for the sin^ 29q determination. 

The error in the measured value of sin^ 29q due to the uncertainty in the value of sin^ 9 is considerably larger 
in an experiment like KamLAND , which is sensitive primarily to the region of the maximum of the De survival 
probability ^ (SPMAX). In this case the oscillatory term s'm^ {AiuqL / iEi,) ^ 1 and can be quite small, so that 
the extra contribution to the uncertainty in sin^20Q, given roughly by Eq. (|30|l . would be rather large. The first 
term in the expression in the right-hand side of Eq. (|30(l becomes relatively large due to presence of the oscillatory 
term in the denominator and the second term becomes even larger due to the presence of the same term both in 
the denominator and in the numerator. Thus, for KamLAND the impact of the uncertainty in sin^ 9 on sin^ 9q 
determination is essentially related to the fact that the Atoq— dependent oscillatory term for KamLAND is relatively 
small, resulting in rather big contributions from both terms in Eq. I|30|l . This explains the relatively large effect of 
the sin^ 9 uncertainty on the sin^ 26*0 determination noticed in • We assert that to reduce the impact of the sin^ 9 
uncertainty on the determination of sin'^6'0, one needs to "tune" the experiment to the SPMIN. This observation 
gives further credence to our statement that the best experimental set-up for determining the solar neutrino mixing 
angle with high precision is a reactor experiment sensitive to the SPMIN. 

In Figs. 121 and 21 we show the regions of allowed values of Amg and sin^6'0, obtained with "data" simulated at 
various values of AttIq > 1.0 x 10"'* eV^ and sin^ 0© — 0.3. We assume a statistics corresponding to an exposure 
of 15 GWkTy (5GWx3kTy) 'si'l, bin size AE = 0.425 MeV and plausible systematic uncertainty of 2%. Figure H 
shows the allowed regions obtained for sin^ 9 — 0, corresponding to the case of 2-neutrino Attiq— driven oscillations. 
In Fig. 01 we display the corresponding allowed regions in a complete 3-neutrino mixing scheme with sin^ 9 allowed 
to vary freely within the 99.73% C.L. allowed range, sin^ 9 < 0.05. As a comparison of the two figures indicates, 
the effect of keeping sin^ 9 free, on the allowed ranges of the solar neutrino oscillation parameters in general, and on 
sin^ 9q in particular, is rather small. This is compatible with the conclusions of the analysis presented above. To 
further substantiate our point pertaining to the impact of the sin^ 9 uncertainty, we show in Fig. a blow-up of the 
99.73% C.L. contours for the sin^ 9 free (within the bound sin^ 9 < 0.05) and the sin^ 9 — cases, obtained with 
L — 30 km, Eth — 2.6 MeV and with assumed "true" value of Atoq — 1.5 x 10~* eV^. We stress that the effect of 
keeping sin^ 9 free is small and henceforth in the rest of this Section neglect the effect of sin^ 9 uncertainty on the 
precision of the Am© and sin^ determination. 

In Fig. 121 we compare the precision obtained on the solar neutrino oscillation parameters for two intermediate 
baselines of 20 km and 30 km, with and without imposing the low energy cut-off of 2.6 MeV ^. For the baseline of 
L = 30 km, depending on the "true" values of the parameters, at 99% C.L. a precision respectively of (6 — 7)% and 
(3 — 6)% is possible to achieve for sin^ 9q and Amg even with Etji = 2.6 MeV. These errors reduce only slightly 



This is valid both for Am| ~ 7.2 X 10"^ eV^ and for Am| ~ 1.5 
* For all the cases we get two solutions in sin^ 6q , the real one and a 
solutions are ruled out by the solar neutrino data. Nonetheless we 



X 10"'' eV^. 

"fake" one on the "dark-side" {sirP 9q > 0.5). This fake "dark-side" 
show them in our plots for the sake of completeness. 




FIG. 3: Simulated two parameter fit of the e"^ —energy spectrum in the ArriQ — sin'^ Oq plane. We assume sin'^ 6* = 0, sin^ 6q = 0.3 
and vary the value of ArriQ. We compare results at 1/ = 20 km and L = 30 km, with the lower cut-off on the energy spectrum 
at Eth ~ 2.6 MeV and at Eth = 1.0 MeV. The energy bin size is chosen as AE = 0.425 MeV, systematic uncertainties are 
taken as 2% and the statistics correspond to 15 GWkTy. Shown are the 90%, 95%, 99% and 99.73% {3a) C.L. contours. 
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FIG. 4: The same as in Fig. |3 but for sin^ 6 allowed to vary freely within its 99.73% C.L. allowed range, sin^ 6 < 0.05. 
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FIG. 5; Comparison of the allowed regions obtained at the 99.73% C.L. for sin^ 9 — and allowing sin^ 9 to vary freely within 
its 99.73% C.L. allowed range, sin^ 9 < 0.05. The curves correspond to the top left-hand panels of Figs. |21 and |1] and "true" 
Ami = 1.5 X 10"" eV^ 



to (5 — 6)% and ~ 2%, respectively, when the threshold is lowered to 1.0 MeV, owing mainly to the fact that the 
statistics increases, but also to the fact that the SPMIN is fully used in the measurement. 

For the L = 20 km experiment on the other hand, the precision for both Atoq and sin^ ^0 (and especially for 
sin^ ^0) is relatively low if the effective threshold of Eth = 2.6 MeV is applied, and if AuIq — 1.5 x 10"'' eV^. As 
discussed in the previous Section, for larger values of Atoq which produce the SPMIN at E^is > 2.6 MeV, it is still 
possible to reach high precision in the determination of sin ^ 6p) and Amg even with the cut-off of Eth = 2.6 MeV. 
These values are somewhat disfavored by the current data [lil [12; HjI • If the threshold is lowered to 1.0 MeV, the 
experiment can "see" the full SPMIN even for Amg = 1.5 x 10"^ eV^ and the sensitivity to both AttIq and sin^ 9q 
becomes remarkable. The AttIq can be determined with a ~ 2% precision, while the sin^ 9q would be known to within 
~ (3 - 4)%, at the 99% C.L. With Eth = 1-0 MeV even though both L = 20 km and L = 30 km baseline experiments 
would be sensitive to SPMIN, the precision of the shorter baseline experiment is higher, owing to the larger statistics 
expected. 

In Fig. 1^1 we study the effect of using smaller bins in the e+— energy spectrum. As explained in the previous Section, 
there is no substantial improvement in the precision of the fit to the solar neutrino oscillation parameters Attiq and 
sin^ 9q with this strategy. The binning energy used by KamLAND {AE — 0.425 MeV) is already sufficient and one 
does not need to use smaller bins. However, we will see in the next sections that a better energy resolution could 
allow for the extraction of other information embedded in the sub-dominant oscillations. 

One of the crucial issues in these types of experiments is thought to be the control of the systematic error. Indeed, 
a large systematic uncertainty on the energy spectrum could threaten to wash out the possibility of extracting any 
information from the data. In Fig. 13 we show the impact of the systematic uncertainty on the expected precision. 
The plots show the allowed areas obtained assuming 1% and 3% systematic uncertainties and can be compared with 
the lower panels in Fig. |21 where the systematic uncertainty is taken as 2%. By comparing the relative precision on 
the parameters, we note that while the precision on Attiq remains essentially unaffected by the exact value of the 
systematic uncertainty if the latter does not exceed ~ 4%, the precision of sin^ ^0 may show a very mild dependence, 
the error in sin^ 9q obviously decreasing with the reduction in the systematic error. However, the impact is not large 
(see also |29j where the performance of KamLAND is studied under assumption of various anticipated systematic 
errors). This implies that, at least for reasonably small values of the systematic uncertainties, their impact on the 
precision of the parameter determination is only marginal. This result is essentially due to the fact that we are 
assuming rather large statistics, allowing for a good reconstruction of the oscillatory pattern in the energy spectrum. 
Similar conclusions were recently drawn in the context of a short baseline experiment at nuclear reactors (L = 1.7 
km) in 

By contrast, we illustrate in Fig.|Slthe huge effect that an increase in the statistical sample could lead to. We show 
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FIG. 6: The same as the lower panels in Fig.|3 but showing the effect of reducing the e^— energy bin size AE. 
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FIG. 7: The same as the lower panels in Fig.|3|but showing the effect of varying the systematic error. 



the allowed regions, corresponding to Fig.O that would be obtained if the statistics were increased 5 fold. This could 
be done using a reactor complex of total power comparable to the Kashiwazaki site (24.3 GW) and a large enough 
detector. We note that the error in the statistical determination of both Amg and sin^ 9q would go down below the 
percent level. With this statistics even the effect of the systematic uncertainty on the precision determination of the 
parameters is expected to be low and one can achieve remarkable accuracy. 
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FIG. 8: The same as in Fig. |3 but showing the effect of increasing the statistics five times, corresponding to 75 GWkTy. 



V. PRECISION PHYSICS FROM Am^^^ -DRIVEN (SUB-DOMINANT) OSCILLATIONS 

In this Section we will consider the set-up where the experiment could be sensitive to the Arn^^.^— driven sub- 
dominant oscillations. As discussed in Section IIIII one has to use relatively small energy bins in order to observe 
these oscillation effects. We first consider the potential of the intermediate baseline reactor experiment of interest in 
improving the bound on sin^ 9. We next assume that sin^ 9 is sufficiently large and can therefore lead to observable 
Am^j,-,-^— driven oscillations and investigate under what conditions one can use these "fast" oscillations to determine 
Am^j,-,-^ with a high precision. Finally, we study the possibility of achieving what is probably a most ambitious goal - 
to get information on the neutrino mass hierarchy by observing the reactor i?e oscillations at intermediate baselines. 



A. Improving the Limit on sin^ 6 

We investigate in what follows the impact of several experimental conditions on the achievable sensitivity on sin^ 9. 
We show in Fig. |^ the sin^ 9 sensitivity obtainable in a 15 GWkTy statistics experiment. We generate the "data" at 
sin^ 9 = and at each value of sin^ 9 we plot the difference between the function obtained at that value of sin^ 9 
and the Xmim which obviously comes at sin^ 9 = 0. We define the 3cr and 90% C.L. limit levels as Ax^ = 9 and 2.71, 
respectively, corresponding to a one parameter fit. 

The sensitivity to sin^ 9 increases substantially if the baseline L is reduced. It follows from Eq. 125|l that the 
best sensitivity to sin^ 9 is achieved for L = few km 113 ■ The intermediate baselines we are considering are not 
optimized for the measurement of 9. The sensitivity to sin^ 9 of the experiment under discussion is worse than that of 
the experiment proposed recently in |4lLl42| with a distance of 1.7 km. Still, the experiment with L = 20 km, would 
allow to put an upper bound of 

sin^ 61 < 0.021 (0.012) , (32) 

at the 3(7 (90%) C.L. (Fig. El top left panel), which is a noticeable improvement over the current 3ct bound of 
sin^ 9 < 0.05. 

Rather relevant for the sensitivity to sin^ 9, in addition to the distance L, is the width of the final state e+ energy 
bins. As could be expected, the smaller bins give a better sensitivity to the value of sin^ 6*. Bins with width of 0.1 
MeV give a 50% improvement in sensitivity compared to the case of bin width of 0.425 MeV. Another relevant factor 
is the presence of a low energy cut-off in the part of the spectrum used in the fit: the sensitivity improves with the 
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FIG. 9: Study of the sensitivity to sin^ 9. Shown are the Ax^ obtained as a function of sin^ 9, when fitting a simulated "data" 
set generated for sin^ 6 = 0, Amatm ~ 2.5 x 10~^ eV^, Amatm = 1-5 x 10~^ eV^ and s'm^ 9q = 0.3. For a particular value 
of sin^ 0, the other parameters are allowed to vary freely in the fit. The two horizontal fines correspond to 3a and 90% C.L., 
respectively. We assume a total exposure of £ = 15 GWkTy, and explore the impact of the baseline (L), width of energy bins 
(AE), the systematic uncertainties and the application of a low-energy cut-off (Eth) on the sin'^ 6 sensitivity. 

inclusion of the lower energy data by, e.g., changing the cut-ofF energy from 2.6 MeV to 1.0 MeV. Finally, the last 
panel in Fig. illustrates the effect of the systematic uncertainties on the sensitivity to sin^6'. As discussed in detail 
in the previous Section, the systematic error in the sin^ 9 measurement has a small but non-negligible effect - the 
limit improving marginally with the reduction of the error. 

We show in Fig. ^| the sin^ 9 sensitivity when the statistics are increased by a factor of five: we consider statistics 
corresponding to 75 GWkTy. Clearly, the increase in statistics improves the s'm^ 9 sensitivity and at 3a (90% C.L.) 
limit for the L = 20 baseline case reads (Fig. El top left panel) 

sin^e < 0.01 (0.0055). (33) 

Thus, for these very high statistics we get a sensitivity to sin^ 9 which is almost of the same order as that expected 
to be reached in the "Reactor-I" experiment discussed in |^. Note that with the increase of the statistics, the effect 
of both the systematic error and of the value of Eth chosen on the sin^ 9 sensitivity decreases. The decreasing of the 
energy bin size, however, has essentially the same effect of increasing the sensitivity. 

B. Measuring Amatm and sin^ 

The experimental setup we consider is optimized for extremely precise measurement of 9q, Atoq, and for distin- 
guishing which hierarchy is realized in the neutrino mass spectrum (as we will see in the next Section), provided the 
high-LMA solution is the correct one. Nevertheless, in order to achieve these goals, a relatively large statistics and a 
sufficiently good energy resolution are required. With these favorable conditions fulfilled and bin width of 0.1 MeV, 
one could hope to improve on the precision of Am^^m ^^'^ ^^^^ provided sin^ 9 is non-zero and is sufficiently large. 
We will show that both of these measurements are rather independent of the value of Arrig, so that they could be 
performed in this type of an experiment even if the low-LMA solution is confirmed by the KamLAND collaboration. 

In Fig. we show the allowed regions in the Am^^jj^ — sin^ 9 plane, obtained for L = 20 km and two different 
statistics samples: the top panels correspond to 15 GWkTy, while the lower panels are for a 75 GWkTy exposure. 
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FIG. 10: The same as in Fig. El but with statistics increased by a factor of five {C = 75 GWkTy). 




FIG. 11: The 90%, 95%, 99% and 99.73% C.L. allowed regions obtained by fitting "data" set simulated at Am^t^ = 2.5 x 10"^ 
eV^ and three different "true" values of sin^ 9. The solar parameters Atuq and sin^ 9q axe taken to be 1.5 x 10~* eV^ and 0.3, 
respectively. The statistics assumed corresponds to £ = 15 GWkTy for the upper panels and to £ = 75 GWkTy for the lower 
panels. The distance used is L = 20 km. The bin size is taken as 0.1 MeV, the lower energy cut-off is Eth ~ 1.0 MeV, and 
the systematic uncertainty is assumed to be 2%. Both the "data" set and the fitted e^— spectrum correspond to NH neutrino 
mass spectrum. 
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FIG. 12: The same as in Fig. 1111 but assuming that the "true" value of Am. 
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The results shown in Fig. are for three different "true" values of sin 9 and for Aml^^-^ = 2.5 x 10^^ eV^, assumins; 
this to be the "true" value. Similar results for Aml^^ = 2.0 x 10~^ eV^ and the same three values of sin^ 9 are 
presented in Fig. El We take a systematic uncertainty of 2% and consider 0.1 MeV bins spanning the entire visible 
energy spectrum between (1.0 — 7.2) MeV. In these figures we keep Auiq and sin^ 9q fixed at 1.5 x lO^"' eV'^ and 0.3 
respectively. Since, as we have shown in the previous Section, the intermediate baseline experiment itself will restrict 
the allowed values of Attiq and sin^ to within a few percent of their "true" value, our results remain practically 
unchanged even if these parameters were allowed to vary freely in the analysis. As Fig.lllldemonstrates. in the 15 (75) 
GWkTy case the experiment under discussion has almost no sensitivity to the value of Am^tm if sin^ 9 ^ 0.02 (0.01). 
Note the similarity between these limiting values of sin^ 9 and the limit on sin^ 9 which the 15 (75) GWkTy experiment 
can provide if the "true" value of sin^ 9 is 0, given in the previous sub-section. This implies that below the indicated 
respective limiting values of sin^ 9, the experiment we are discussing cannot distinguish between a zero and a non- 
zero sin^ and thus fails to "see" the Arn^t^— driven oscillations. The reason for the modest sensitivity to sin^ is 



'atm 

"our' 

oscillations. Correspondingly, L is much bigger than the L* (cf. Eq. 1)25(1 ') 



related to the fact that the L is chosen for "our" experiment to be sensitive to the SPMIN due to the AuIq 



driven 

suitable for the SPMIN associated with the 



Am^j,-,-^— driven oscillations. The lack of accuracy in the reconstruction of the details of a single oscillation period (in 
Evis) affects strongly the sensitivity to sin^ 9 which controls the amplitude of the oscillations. However, this should not 
affect the sensitivity to Am^^j^j, which determines the oscillation length, because of the presence of many oscillation 
periods in the relatively large energy window offered by the spectrum itself. 

Thus, if sin^ 9 is large enough so that the Ato^^j^— driven oscillations can be observed in the experiment under 
discussion, Am^^j^^ can be measured with a high precision. Indeed, in the case when the 
Ami 

the sensitivity of the JPARC (JHF-SK) project to Arn^ 

^atm ' 

the L = 1.7 km experiment is optimized to see the Am^^j^— driven oscillations. However, one should keep in mind 
that the L — 1.7 km baseline experiment cannot constrain the solar neutrino oscillation parameters. The uncertainty 
in ArriQ allows Atoq to take on values as high as (2 — 3) x lO^** eV^. This drastically reduces the sensitivity of the 



'true" value of sin 9 = 0.03, 

can be determined to a within few percent accuracy at the 99% C.L. This accuracy could be comparable to 

project to Am'L^ ^ It is certainly much better than the sensitivity to 
Amitry, of the L = 1.7 km short baseline experiments [ly with the same statistics. This is somewhat unexpected since 



^ The experiment we are discussing has a much larger statistics than the planned JPARC (JHF-SK) experiment; for the precise determi- 
nation of Am^jj^ it also requires sin^ 9 to be sufficiently larg 
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For this uncertainty the authors of _4,lJ have used the current 3cr uncertainty on Am; 
KamLAND data analysis. 
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FIG. 13: The impact of the "true" value of Attiq on the 99.73% allowed areas when the "true" Amatm = 2.5 x 10 eV^ and 
sin^ 9 = 0.03. 

L = 1.7 km experiment to Am^^j^^. In the experimental set-up we are discussing, the solar parameters are determined 
to a very high degree of accuracy by the experiment itself and no external input on their errors is required. Thus, 
the intermediate baseline experiment has a better sensitivity to Am^^j^, provided sin^ 9 is sufficiently large, while the 
sensitivity to sin^ 6 of the shorter baseline experiment is better for the reasons discussed above, viz. it is optimized 
to the Am^^j^-driven SPMIN. 

Finally, we show in Fig. ^| the impact of the "true" value of Atoq on the ability of the experiment under con- 
sideration to measure Arn^^^ and sin^ 9. The figure very clearly demonstrates that even if the low-LMA solution is 
confirmed as the true solution of the solar neutrino problem, the intermediate baseline experiment could still be used 
to measure Am^j^ and/or constrain sin^ 6. 

We stress that Figs. ITTl [HI and [HI are obtained using Eq. JHl, i.e., assuming that the neutrino mass spectrum 
conforms to a NH. Our conclusions for the sensitivity to sin^ 9 and Aml^^ would remain the same for the case of 
IH mass spectrum. For a given "data" set, the allowed region in the parameter space, especially the fitted values of 
Am^jjjj, depend on the hierarchy assumed if the Attiq is in the high-LMA region. The possibility of two types of 
neutrino mass hierarchy leads to an ambiguity in the allowed values of Am^^j^^ obtained, as we will show in the next 
sub-section. However, even though the assumption of the "wrong" hierarchy can lead to another (separated) allowed 
region in the parameter space, it does not describe the data so well as the "correct" hierarchy and is disfavored, in 
general. This can be used to gain insight into the type of hierarchy the neutrino mass spectrum has and can allow 
"our" intermediate baseline experiment to determine the hierarchy. 

Let us stress that the question of the hierarchy becomes relevant if Amg lies in the high-LMA zone. This will be 
the subject of our discussion in the next subsection. If the true solution turns out to be low-LMA, the last terms in 
Eq. jnj and Eq. (|l()|l would be negligible and sensitivity to hierarchy would be lost. Let us emphasize, however, that 
even in this case it would still be possible to improve the existing limits on sin^ 9, and - if sin^ 9 is sufficiently large - 
to measure Arri^^^ with an exceptional precision. 

C. Normal vs. Inverted Hierarchy 

The next question that may be answered with the experiment under discussion is whether the neutrino mass 
spectrum is with normal hierarchy or with inverted hierarchy. As pointed out in Section fllll and follows from Eq. H21|l . 
the baseline L = 20 30 km is particularly suited for that purpose. In this subsection we present results for L = 20 
km. We use the highest statistics and smallest energy bins and consider the entire visible energy spectrum with a 




FIG. 14: The energy spectrum (in 0.1 MeV bins) of the events in the case of C-e oscillations (left-panel), and the ratio of events 
in the cases of oscillations and absence of oscillations (right-panel), showing the effect of hierarchy for a L = 20 km experiment. 
The thick cyan/grey line corresponds to the case for NH with Arriatm = 2.5 x 10"'' eV'^ while the dotted and thin solid line 
correspond to the case of IH with Am^tm ~ 2.5 x 10~^ eV^ and Amatm = 2.6 x 10"'^ eV^, respectively. The statistics assumed 
is £ = 75 GWkTy, while sin^ 9 = 0.03, Am% = 1.5 x lO"" eV^ and sin^ 9e = 0.3. 



low energy cut-off of Eth ~ 1.0 MeV. We assume also that the "true" value of sin^ 6 is non-zero and is sufficiently 
large. As it follows from Eq. © and Eq. 1)10(1 . the predictions for the final state e^— spectrum distortions in the 
NH and IH cases differ only if the solar neutrino mixing angle Oq ^ 7r/4. Maximal mixing is currently disfavored by 
the solar neutrino and KamLAND data [Til . We use the current best-fit value for the solar neutrino mixing 

angle, corresponding to sin^ 9q = 0.3. The value of sin^ 6q will be measured with a very high precision in the reactor 
experiment under discussion itself. 

In Fig. ^] we show the visible energy spectrum expected for both the NH and IH cases. The left-hand panel shows 
the number of events, while the right-hand panel shows the ratio of the number of events in the case of De oscillation 
("oscillation events") to the number of events in the absence of oscillations ("no oscillation events"). The thick 
cyan/grey line is for the case of Aml^^ = 2.5 x 10~'^ eV^ and NH, the dotted fine corresponds to Am^^j^^ = 2.5 x 10~^ 
eV^ and IH, while the thin solid line is for Am^tm = 2.6 x 10""^ eV^ and IH. Thus, for the same value of Am^^^, the 
6+— spectrum deformations expected for NH and IH are different owing to the last terms in Eq. © and Eq. 1(10(1 . One 
might expect on the basis of Eqs. © and ((10(1 that the IH e"*"— spectrum would fit a NH e+— spectrum with a value 
of Am^jjjj that is larger than the "true" value by a range which would be of the order of Ami,. If the experiment 
has a sensitivity to Am^^jjj which is better than, or is at least of the same order as Attiq, one might expect two 
non-degenerate solutions in the Arnl^^ — sin^ 6 parameter space for the same data set: one for NH and another for 
IH. 

However, we would like to stress that even though the IH spectrum could approximately reproduce the NH 

e"^ — spectrum, it cannot exactly reproduce the latter. Hence, if the NH is the true hierarchy, the case of IH mass 
spectrum would always be disfavored by the data from "our" experiment compared to the NH spectrum. We have seen 
in the previous Section that the experimental set-up we are discussing could have a very high sensitivity to Am^jj^j, 
and therefore we can expect the large statistics intermediate baseline reactor experiments to determine the type of 
neutrino mass hierarchy, for sufficiently large sin^ 9 at least, if the high-LMA solution holds. 

In Fig. El we explicitly show the allowed regions obtained using a "data" set generated in the NH case and fitted 
by both the NH and IH expressions for the Vg survival probability. The results shown are obtained using the 75 



0.003 
0.0029 
0.0028 
0.0027 



0.0026 
^ 0.0025 
I 0.0024 

0.0023 

0.0022 

0.0021 
0.002 

0.01 0.02 0.03 0.01 0.02 0.03 0.04 
sin^S sin"6 



FIG. 15: The same as in Fie. II II but probing the type of the neutrino mass hierarchy. The "data" sets simulated correspond to 
NH. The filled contours show the C.L. allowed regions obtained by fitting for the parameters assuming NH. The empty contours 
show the 90%, 95%, 99% and 99.73% solution regions, obtained by fitting the NH "data" with an IH theory. The contours for 
the IH are obtained with respect to the x'^-minima found in the case of the NH. The statistics assumed in the left-hand panel 
is £ = 75 GWkTy, while that in the right-hand panel is £ = 125 GWkTy. We note that while for the £ = 75 GWkTy case the 
"wrong" hierarchy is allowed even at the 90% C.L., with the L = 125 GWkTy statistics the experiment can rule out the IH at 
least at the 95% C.L. 

GWkTy statistics and a higher 125 GWkTy statistics. The "true" value of Am^^j^j assumed is 2.5 x 10"'^ eV^, while 
"true" sin^ 6* = 0.02. For the 75 GWkTy case, the two regions overlap at the 3tT level, while for 125 GWkTy they 
are completely non-degenerate. Note also that for the 125 GWkTy case, the "data" generated for NH spectrum can 
"rule out" the IH spectrum, at least at the 95% C.L. 

In Fig. 115! we show the dependence of the hierarchy sensitivity of the L = 20 km reactor experiment on the "true" 
values of Am^^j^, sin^ Q and Atoq. For each point in the Amatm — sin^ Q parameter space and a given Amg and 
sin^ 00, we simulate the observed e"'"— energy spectrum assuming NH to be true. We then fit this "observed" spectrum 
with both the NH and IH through a analysis, allowing the parameters Ato^^jj^ and sin^ Q to vary freely around 
their "true" values. The IH spectrum obviously does not fit the "data" as well as the NH and hence is disfavored. To 
quantify this sensitivity to hierarchy we calculate the difference (x^iff) between the for the NH and IH for each 
point in the parameter space and use the definition Xdiff < to find whether IH can be ruled out at a given C.L. 
The Ax^ is determined by the C.L. considered and we take the 2-parameter definition for it. We show the lines of 
the C.L. contours in the Arnl^^ — sin^ 9 plane for 2 different values of Ami, and sin^ 6q = 0.3. All "true" values of 
the parameters located to the right of a given line, could allow the experiment to disfavor the IH spectrum at the 
corresponding C.L. The solid lines give the sensitivities for 75 GWkTy exposure, while the dashed lines are for a 125 
GWkTy data set. The sensitivity to hierarchy is maximal for the smallest allowed values of Am^^j^ and for these 
cases it is possible to rule out IH even for smaller values of sin^ 9. It becomes increasingly more difficult to distinguish 
between the two types of neutrino mass hierarchy as Arn^^^^ increases, as the figure indicates. The sensitivity of the 
experiment to hierarchy depends critically on the "true" value of the solar parameters, as is evident from Eq. @ and 
Eq. (|10|l . In particular, higher values of Atoq are better suited for hierarchy determination, as can be seen explicitly 
comparing the two panels in Fig. 116! 

The possibility of determining the neutrino mass hierarchy is very sensitive to the value of sin^ 9q . As we have 
explained above, when we generate the NH "data" and fit them with the IH theory, we allow Am^j,-,-^ to vary and take 
larger values than Ato^^j^ used to generate the "data". Let us denote the best-fit Ato^^j-^ in the IH case as Am'^^^^-^: 
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FIG. 16: Contours of constant 99.73%(black), 99%(red), 95%(green) and 90%(blue) C.L., showing the abihty of the L = 20 km 
baseline experiment to determine the correct hierarchy at each "true" point in the Am^tm — sin^ 9 plane. The solid lines (dashed 
lines) correspond to £ = 75 GWkTy (£ = 125 GWkTy). The resuhs shown are for sin^ Oq = 0.30, and for Am|, = 1.5 x 10"* 
eV^ (right panel) and Ajtiq = 2.5 x 10~* eV^ (left panel). For the values of Am^t^j and sin^ 6, corresponding to points on the 
plane, located to the right of a given C.L. line, the experiment can rule out the "wrong" hierarchy at that C.L. See the text 
for further details. 



where Am > 0. The difference between the IH and NH survival probabilities at each E is given by: 



PlH(Ve — > Ve) - PNH{Ve ^ Vf.) = 4 sin^ 6* COS^ 



' cos^ 6*0 sin( ^^2^— + ^) sin( ^ 

° ^ 2E 2E AE ' ^ AE ' 

-sm^^esm(^f^-^)sn.(^) 



(35) 



When, e.g., cos 26*0 — 1, we have Pjh — Pnh — for Am? = AttIq. For cos29q ~ (maximal mixing), the difference 
between the two probabihties is for Am? = 0. For the reahstic case of cos2^^0 = (0.10 — 0.40), it is impossible to 
have PiH — -Pnh = for every value of E from the interval of interest. This is illustrated in Fig. El where we show the 
difference between the minimal values of for the "wrong" IH and the "right" NH spGctra- cis a function of sin 6q 
for "true" Am| = 1.5 x lO'^; 2.5 x 10"'*; 3.5 x 10"^ eV^, Am^^,^^ = 2.5 x lO'^ eV^ and sin^ 6* = 0.03. As Fig. El 
indicates, the best sensitivity to the type of neutrino mass hierarchy is achieved in the interval sin^ Oq = (0.3 — 0.45), 
which is favored by the current solar neutrino and KamLAND data. The rather strong dependence of the sensitivity 
to the type of neutrino mass hierarchy on sin^ 9q is illustrated in Fig. 18, in which we show contours of constant C.L. 
in the Am^^^-^ — sir? 9 plane, at which one could distinguish between the NH and IH spectra, for Attiq = 3.5 x 10"'' 
eV2 and sin^ Oq = 0.30; 0.40. 

One final comment is in order. As Fig. 1151 illustrates, for a given "data set" (i.e., "true" value of Am^^^-^, etc.) and 
sufhcient statistics, the neutrino mass hierarchy ambiguity can lead to two non-degenerate solutions for Am^^j^^, one 
for NH and another for IH spectrum. In certain of these cases the intermediate baseline reactor experiment might not 
be able to rule out the "wrong" hierarchy even at the 90% C.L. (Fig. El left panel). However, since the NH and IH 
give two different and non-degenerate values for Am^^.^, in such cases one can use the information from some other 




FIG. 17: The difference between the minimal values of x^, obtained by fitting a NH "data" set with NH and IH spectrum, 
as a function of sin^ 9q . The results shown are for Am| = 1.5 x 10"^ 2.5 x 10"^ 3.5 x 10~^ eV^ The "true" values of 
Am^tm = 2.5 X 10"^ eV^ and of sin^ 9 = 0.03, while the statistics corresponds to £ = 75 GWkTy. 




FIG. 18: The same as in Fig. [TSl for £ = 75 GWkTy, Am| = 3.5 x 10"* eV^ and two values of sin^ Gq: 0.30 (solid lines) and 
0.40 (dash-dotted fines). 



very precise experiment, like the proposed JHF-SK or NuMI ofF-axis experiments with neutrino superbeams |43l l44| . 
to determine the correct value of Am^^^. The latter can be used together with the data from the reactor experiment 
to get information on the hierarchy. 



VI. CONCLUSIONS 



In this paper we analyzed the physics potential of a reactor neutrino experiment with a relatively large detector 
at a distance L ~ 20 30 km. This distance has been chosen in order to achieve the best sensitivity to the solar 
neutrino oscillation parameters Amg and sin^ 9q , assuming that the latter lie in the high-LMA solution region. We 
considered the case of three flavor neutrino mixing and used in the analysis the exact expression for the relevant 
j>e survival probability. The latter depends, in addition to Auiq and sin^^©, on Am^^j^^, driving the atmospheric 
neutrino oscillations, on 9 - the angle limited by the CHOOZ and Palo Verde experiments, and for sin^ 9q < 0.50 - 
on the type of the neutrino mass spectrum which can be with normal hierarchy (NH) or with inverted hierarchy (IH) . 
The current solar neutrino and KamLAND data favor sin^ 9q ^ 0.30. We discussed strategies and the experimental 
set-up, which would permit to measure Atoq and sin^ 9q with a high precision, get information on (or even measure) 
sin^ 0, and if sin^ 9 is sufficiently large (sin^ 9 <; 0.02) provide a high precision measurement of Am^tm and determine 
the type of the neutrino mass hierarchy. More specifically, we have investigated the impact that i) the choice of the 
baseline L, ii) the effect of using a relatively low e"*"— energy cut-off of Eth ~ 1-0 MeV, iii) the detector's energy 
resolution, as well as iv) the statistical and systematical errors, can have on the measurement of each of the indicated 
neutrino oscillation parameters. 

The precision with which Atoq and sin^ 9q can be determined in the experiment under discussion depends crucially 
on whether the first Amg— driven oscillation minimum of the Ue survival probability falls or not in the interval of 
e"*"— energies which can be used for the measurements in the experiment. For Amg = 1.5 x 10^"* eV^ (2.5 x 10^** eV^) 
and L = 20 km, this minimum takes place at the e"*"— energy of Eyis = 1.6 (3.2) MeV; for L = 30 km, it occurs 
at Eyis = 2.8 (5.3) MeV. This implies that if Am© = 1.5 x 10~^ eV^, in order to achieve a high precision in the 
determination of Atoq and especially of sin^ 9q in an experiment with a baseline of L = 20 km, a relatively low 
"threshold" energy should be employed, Eth ~ 1-0 MeV. If, however, Atoq ^ 2.0 x 10~* eV^, one can use Eth — 2.6 
MeV, as was done in the KamLAND experiment. In the case of L = 30 km one can use Eth — 2.6 MeV even for 
Am© ^ 1.5 X 10"'' eV^. If the condition under discussion is fulfilled, remarkable precision in the measurement of 
Am© and sin^ 9 can be achieved. For the lower "threshold" energy, Eth ~ 1-0 MeV, achieving high precision requires 
a controllable or negligible background due to the geophysical neutrinos in the interval Eyis = (1.0 — 2.6) MeV. If 
a sufficiently accurate modeling of this background is achieved and Eth — 1-0 MeV, assuming Atoq = 1.5 x 10"'' 
eV^, energy bins of AE = 0.425 MeV, sin^ 9q = 0.3, 2% systematic error and a total statistics corresponding to 
£ = 15 GWkTy, one could determine Atoq with up to a ^ 2%, and sin^ 9q with a 3 4% uncertainty at 99% C.L. 
in an experiment with a baseline of L = 20 km (Figs. 5). This impressive results depend rather mildly on the 
systematic error provided the latter does not exceed ~ 4% (Fig. Ej). The reduction of the e"*"— energy bin width leads 
to a negligible change in the precision with which Atoq and sin^ 9q are measured (Fig. O . The precision can be even 
higher in the case of a larger statistics (e.g., for £ = 75 GWkTy, Fig.lSj). 

The same apparatus could be used to extract information on the other parameters entering into the expression for 
the survival probability, through the study of sub-leading effects controlled by the (small) mixing angle 9, for which 
only upper bounds exist at present. 

The limit on sin^ 9 can be lowered to sin^ 9 < 0.021 at 99.73% C.L. for L = 20 km, if the "threshold" energy 
Eth — 1-0 MeV, the systematic error is 2%, the e"*"— energy bin width is AE = 0.1 MeV and if statistics corresponding 
to £ = 15 GWkTy is collected. This result depends strongly on the choice of L: the optimal distance for such a 
measurement is of a few km, and the shorter the distance with respect to L = 20 km the stronger this bound would 
be. The magnitude of e"*"— energy bin width and the value of the e"*"— energy effective threshold have a considerable 
impact on the bound as well, while a further reduction of the systematic errors could lead only to a mild improvement. 
If the statistics is as large as that corresponding to a set up with £ = 75 GWkTy, one could improve the upper bound 
to sin^ 9 < 0.010 (0.0055) at 99.73% (90%) C.L. 

If, on the contrary, a relatively large non- vanishing value of sin^ 9 is found, the experiment could gain sensitivity to 
Am^tm through the detection of the sub-leading oscillatory properties of the survival probability. Not averaging out 
the Arn^jjjj— driven oscillations is crucial for the measurement of Arn^^.^. This requires a relatively high e"''— energy 
resolution, permitting a binning in the e'*'— energy with a rather small bin width, AE = 0.1 MeV. Taking Eth = f -0 
MeV and assuming 2% systematic error, we find that a set-up with £ = 15 (75) GWkTy could provide a precise 
determination of Aml^^ with a few percent uncertainty for Aml^^ 2.5 x 10~^ eV^, if sin^ 9 ~ 0.03 (0.02) (Fig. [TT|l . 

Finally, the distance L = 20 km is the optimal one for distinguishing between the NH and IH neutrino mass 



spectrum. The relevant v,, survival probabilities for the two types of spectrum differ for sin^ Oq < 0.5 by a 
sin^ 0— suppressed interference term due to the amplitudes of the Amg— and Arn^^^— driven oscillations. For this 
highly challenging study, a very large statistics is necessary: we considered £ = 75 and 125 GWkTy. The possibility 
of distinguishing between the two types of neutrino mass hierarchy was found to depend not only on the values of 
sin^ 9 and sin^ 9q, but also on the values of Amg and Am^^j^^. The dependence on Attiq and sin^ 6q is particularly 
strong. For Amg = 1.5 x 10"'' eV^, and C ~ 75 (125) GWkTy, one could distinguish the NH from the IH spectrum at 
99.73% C.L. in the region of Aml^^, < 2.5 x lO'^ eV^ if sin^ 9 > 0.038 (0.03) (Fig.Cni). The type of the neutrino mass 
spectrum can be determined for considerably smaller values of sin^ 9 if Atoq has a larger value (Figs. 17 and 18), or 
if an even larger statistics would be accumulated. For £ = 75 GWkTy and Atoq = 2.5 x 10~^ eV^, for instance, this 
can be done at 99.73% C.L. for sin^ 9 > 0.019 if Am^^^ < 2.5 x 10"^ eV^ and sin^ 9q = 0.30 (Fig. As we have 
shown, the highest sensitivity to the type of the neutrino mass hierarchy is achieved for sin^ 9q = (0.30 — 0.45) (Figs. 
17 and 18). 

To conclude, the intermediate baseline L = 20 km reactor neutrino experiments have in the case of the high-LMA 
solution of the solar neutrino problem a remarkable potential for a high precision measurement of the solar neutrino 
oscillation parameters Amg and sin^ 9q, which can reduce the error in the values of the latter to a (2 — 4)%. Such 
an experiment can also improve the existing bounds on, or measure, sin^ 9, currently limited by the CHOOZ and 
Palo Verde data. If sin^ 9 turns out to be relatively large, sin^ 9 ^ 0.02, the same experiment could measure Aml^^^, 
driving the atmospheric neutrino oscillations, with an error of a few percent and might be able to establish whether 
the neutrino mass spectrum is with normal or inverted hierarchy. 
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